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With an inferred bolometric luminosity of 5 x 1O 15 L , the recently identified 
z=3.87, broad absorption line quasar APM 08279+5255 is apparently the most 
luminous object currently known. As half of its prodigious emission occurs 
in the infrared, APM 08279+5255 also represents the most extreme example 
of an Ultraluminous Infrared Galaxy. Here, we present new submillimeter 
observations of this phenomenal object; while indicating that a vast quantity 



(N 



O of dust is present, these data prove to be incompatible with current models of 

emission mechanisms and reprocessing in ultraluminous systems. The influence 
of gravitational lensing upon these models is considered and we find that while 
the emission from the central continuum emitting region may be significantly 
enhanced, lensing induced magnification cannot easily reconcile the models with 
| observations. We conclude that further modeling, including the effects of any 



c3 



differential magnification is required to explain the observed emission from APM 
08279+5255. 
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1. Introduction 

With bolometric luminosities exceeding 1O 12 L , Ultraluminous Infrared Galaxies 
(ULIRGs) represent an extreme class of objects whose spectra are dominated by emission 
in the far-infrared (for a review of their properties see Sanders fc Mirabel 1996|) . The 



source of their prodigious output is thought to arise in a thick, cool (T~100°K) nuclear 
dust structure which reprocesses emission from an obscured AGN-like core, a massive 
star-formation region, or possibly both ( Genzel et al. 1998 ). The tunneling of fuel into 



the nucleus and the resultant activity may be triggered by a violent interaction between 
gas-rich galaxies ( [Barnes fc Hernquist 199"6| , [laniguchi fc Shioya 1998] ) ; this appears to 



be the case observationally, with ~ 90% of ULIRGs displaying disturbed or merging 
morphology ( |Clements fc Baker 1996| ). 



Recently identified in a survey of halo carbon stars, the z=3.87 Broad Absorption Line 
(BAL) quasar, APM 08279+5255, was found to be positionally coincident with a source in 
the IRAS Faint Source Catalog fllrwin et al. 1998|) . With a flux of 0.9Jy at 100/mi, this 
m r = 15.2 object possesses an inferred bolometric luminosity of 5 x 1O 15 L (Q = 1 and 
H = 50km/s/Mpc throughout), making APM 08279+5255 apparently the most luminous 
object currently known. However, ground-based images taken in ~ 0".9 seeing reveal that 
the source is slightly elongated, suggesting that APM 08279+5255 consists of a pair of 
sub-components separated by ~ 0'.'3, consistent with the merger-driven hypothesis but also 
with the action of gravitational lensing. In several extreme ULIRGs the morphology of the 
sub-components directly reveals the action of gravitational lensing, the magnification effect 
of which can significantly enhance the intrinsic properties of a system; this is the case in 
both H1413+117 (the "Cloverleaf Quasar" |Magain et al. 1988| , |Kneib et al. 1998Q and the 



hyperluminous galaxy IRAS F10214+4724 ( [Rowan-Robinson et al. 199 1| , Broadhurst fc 



[Lchar 1995 , Eiscnhardt et al. 1996 ), where the lensing-induced amplification is estimated to 



be ~ 30 — 100. As APM 08279+5255 is apparently so extremely luminous, the possibility 
that gravitational lensing is influencing the observed properties is highly likely. 

Given the cool temperature of the dust region, photometric and spectroscopic 
observations at submillimeter wavelengths have proved to be an important probe of ULIRG 
systems, revealing details of the physical properties and processes underway in these 
energetic objects (e.g. |lvison et al. 1998a] , [Ivison et al. 1998"5| , [Hughes fc Dunlop 1998| ). In 



this paper, we present new submillimeter photometry of APM 08279+5255, from which an 
estimate of the total mass of dust, and its temperature, are derived. Combining these data 
with previous observations, we compare the total spectral energy distribution to popular 
models of ULIRG galaxies, while taking into account the potential effects that gravitational 
lensing may be playing in distorting our view of this system. 
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2. Observations 



The observations were conducted with the Submillimeter Common- User Bolometer 
Array (SCUBA, |Cear fc Cunningham T&JE) on the James Clerk Maxwell Telescope 



SCUBA contains a number of detectors and detector arrays cooled to 0.1°K that cover 
the atmospheric windows from 350/xm to 2000/im. For our photometric observations, we 
operated the 91 element Short-wave (SW) array at 450/mi, the 37 element Long-wave 
(LW) array at 850/zm, and the single photometry pixel at 1350/im, giving half-power beam 
widths of 7.5, 14.7, and 21" respectively. A 9-point jiggle pattern was employed to reduce 
the impact of pointing errors. The source was centered on the central pixel of the arrays 
with the outer pixels being used for subtraction of the sky variations. Whilst jiggling, the 
secondary was chopped at 7.8125 Hz by 60" in azimuth. The pointing stability was checked 
every hour and regular skydips were performed to measure the atmospheric opacity. The 
rms pointing errors were below 2". 

The observations were conducted over 3 nights in April, 1998. The first night had 
stable atmospheric zenith opacities at 450/im and 850/im, with r being 1.47 and 0.29, 
but extremely high winds which forced dome closure. Firm detections were obtained 
at both wavelengths. Sky conditions were reasonably good on the second night and 
observations at both 850/im/450/im(r = 0.58/3.7) and 1350/xm were carried out. On the 
third night conditions were less favorable (CSO r around 0.1 and variable) and only 1350/im 
observations were done Q. 

On the first and second nights both CRL618 and IRC10216 were observed as calibrators, 
while only IRC10216 was observed on the third night. The latter is variable on a timescale 
of two years, and all observations were referenced to CRL618, which has well-determined 
flux densities with SCUBA at 850/zm and 450/zm, but nothing yet published for SCUBA 
at 1350/im. In the latter case, the CRL618 flux density published by Sandell (1994) was 
used, and resulted in a sensible value for IRC10216. The values derived for IRC10216 at 
450/im/850yum from CRL618 were also well within the quoted errors. 

The dedicated SCUBA data reduction software (SURF, |Jennes fc Lightfoot 1998Q was 
used to reduce the observations. All non-noisy bolometers beyond 40" of the central pixel 



6 The James Clerk Maxwell Telescope is operated by The Joint Astronomy Centre on behalf of the Particle 
Physics and Astronomy Research Council of the United Kingdom, the Netherlands Organisation for Scientific 
Research, and the National Research Council of Canada. 

7 Observations on days 2 & 3 were undertaken as part of the CANSERV program operated by the Herzberg 
Institute for Astrophysics. 
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were used to compensate for spatially-correlated sky emission in the 850//m/450/im arrays. 
The 1350/xm pixel currently has no provision for subtracting sky variations using the other 
wavelength pixels. The results of these observations, as well as variance-weighted mean 
values, are listed in Table |I| 



3. Spectral Energy Distribution 

The submillimeter photometry presented in Table [I] is combined with previous data 
(Irwin et al. 1998) in Figure [TJ, which shows the spectral energy distribution (SED) of APM 
08279+5255. The new data reveals that, in the submillimeter regime, the SED possesses 
a slope of 3.1 ± 0.2 in vF u , consistent with Rayleigh- Jeans blackbody emission. The 
curves superimposed on Figure [l] show pure blackbody spectra of a source of temperature 
T=120°K(dot-dash line) and T=220°K(dotted line); such temperatures are representative 
of the temperature range of the regions of the source that emit in the submillimeter /far 
infrared region of the SED. These values are slightly higher than those determined for the 
far-infrared emission in other ULIRG systems, which are typically T ^ 100°K ( [Kales fc 



Edmunds 1996| ). The implied radius of such a blackbody is ~ 650pc (assuming the emission 



region is spherical), a factor of 2 smaller than the cooler (80°K) region responsible for the 
emission in IRAS F10214+4724 ( Pownes et al. 1992Q . 



As outlined in Hildcbrand (1983~| , observations in the submillimeter can be used to 



determine the mass in dust of ULIRGs (e.g. |Clements et al. 1992| , |Downes et al. 1992] , [Eales| 



Bz Edmunds 1996| ). Typically, the dust emission in these systems appears to be optically 
thin, as is inferred from their steep SED slope (the blackbody spectrum is modified by 
a dust absorption coefficient, Kd oc u n , where n = 1 — 2; in the Rayleigh- Jeans regime, 
this is apparent as a steepening of the slope to at least 4 in vF v ). The SED of APM 
08279+5255 does not show such a steep slope, it is instead consistent with emission from an 
optically thick blackbody. The dust mass obtained from the application of the optically thin 
treatment can, therefore, be treated as a lower bound on the total dust mass in the system. 
Assuming that the dust has a temperature of 220°K, the resulting dust mass from the three 
submillimeter data points presented here is 3.7 + 0.3 x 1O 9 M [assuming a constant dust 
absorption coefficient of k& — 0.1 m 2 /kg ( [Hughes fc Dunlop 1998] )]. If gravitational lensing 
plays a role, this value could be over-estimated by a factor of ~ 30, although it does indicate 
that APM 08279+5255 possesses copious amounts of dust, equivalent to that inferred for 
IRAS F10214+4724 ( Pownes et al. 1992| ), and in other ULIRGs ( [Eales fc Edmunds 1996| ). 



The complex form of the SED in ULIRGs suggests that a single temperature 
blackbody source for the FIR flux represents a gross simplification of the underlying 
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processes, and several multi-component models have been developed to better model 
these systems ( [Rowan- Robinson et al. 1993| , [Efstathiou fc Rowan- Robinson 1995| , preen fc| 
Rowan- Robinson 1996| , |Granato et al. 1996| ). These models generally consider a powerful 



source of continuum radiation embedded in an thick distribution of dust possessing grains 
of varying composition ( [Rowan- Robinson 1986|) . Radiative transfer techniques are then 
employed to calculate the reprocessing and re-emission of the central continuum radiation. 

Three such models are superimposed on the submillimeter-optical region of the 
SED in Figure |^. The dot-dashed curve represents the "Embedded Quasar Model" 
of [Rowan- Robinson et al. 1993| (their Model C); this possesses a power-law continuum 



source at the center of a spherically-symmetric dust distribution. As the dust obscures a 
direct view of the illuminating source, all observed emission has been reprocessed. This 
model spectrum is normalized to the submillimeter data points presented here. Such a "pure 
dust emission" spectrum provides a poor fit to APM 08279+5255 and other ultraluminous 
systems, severely underestimating the optical flux. Spherical models with a lower optical 
depth allow some continuum flux through the dust region, but such models cannot account 
for the far infrared/sub millimeter flux ( [Rowan- Robinson et al. 1993 ). 



A natural extension to such a model is to remove the assumption of spherical symmetry 
and consider instead an axisymmetric distribution of absorbing dust, representing a torus 
about the nuclear region ([Rowan-Robinson et al. 1993| , [Efstathiou fc Rowan-Robinson 



1995| , preen fc Rowan- Robinson 1996|) . As well as scaling with total luminosity, the 



emergent spectrum is also then a function of orientation with respect to the observer; 
when viewing a system from the equatorial plane the dusty torus obscures the view of the 
continuum source and the SED is dominated by infrared dust emission. When viewing 
from the pole, however, an observer can look directly onto the central continuum source 
and the SED can possess both significant infrared and optical components. Such a model 
was recently applied to several of extreme ultraluminous systems and, considering only 
orientation and total luminosity, was found to reproduce the gross characteristics of the 
observed SEDs ( Granato et al. 1996|) . The solid line in Figure ^| of this model, as viewed 



from the pole, while the dotted line is an equatorial view. Also plotted are the rest frame 
SEDs for the ultraluminous Cloverleaf quasar, H1413+117 (filled triangles), and IRAS 
FSC 10214+4752 (crosses), both normalized to the submillimeter observations of APM 
08279+5255. While adequately describing these spectra, it is apparent that the models of 
Granato et al. underestimates the infrared to optical flux by a factor of £ 5 even at the 
extreme polar viewing orientation, and provide a poor fit of the SED of APM 08279+5255. 
Extrapolating this model to the data presented in this system, the implied dust mass is 
~ 3 x 1O 8 M , but again likely underestimates the true mass of dust in this system and is 
dependent on the degree of gravitational lensing. 
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4. The Influence of Gravitational Lensing 

The degree to which the flux from a source is enhanced by the action of gravitational 
lensing is dependent on the scale-size of the source, leading to pronounced differential 
magnification effects flSchncidcr et al. 1992 ); in IRAS F10214+4724 the optical emission is 



thought to be magnified by a factor of ~ 100, while the flux from the more extended far- 
infrared region undergoes a magnification of ~ 30 (|Eisenhardt et al. 1996| ). If gravitational 



lensing influences APM 08279+5255, can such effects account for the discrepancy between 
these models and the observed SED? In "typical" quasars, the optical-to-infrared continuum 
ratios are seen (in vF v ) to be ~ 2 flSaunders et al. 1989| , [Elvis et al. 1994] ). All this 



emission arises in the central accretion disk whose scale is similar for both emission regimes; 
gravitational lensing would uniformly magnify such a source, simply scaling its contribution 
to the total spectrum. Normalizing the (magnified) quasar SED to the observed R-band 
magnitude it can be seen that the quasar continuum source, even if highly magnified in 
this system, can contribute little to the submillimeter-infrared SED. It should be noted, 
however, that if the quasar continuum source does, intrinsically, possess significant emission 
into the far-infrared, differential magnification effects may reconcile the data with the 
current models. 



5. Conclusions 

This paper has presented new submillimeter observations of the ultraluminous 
BAL quasar APM 08279+5255; the data are consistent with emission from a warm 
(120°K-220°K), massive, optically thick distribution of dust which is heated by a quasar 
central continuum source. A simple model for the emission region of ULIRGs, consisting of 
a AGN-like continuum embedded within a spherical distribution fails, when normalized to 
the submillimeter data, to reproduce the observed SED, underestimating the the flux in 
the infrared and, due to a highly obscured view onto the quasar source, predicts no optical 
flux. Axisymmetric models, where the dust is distributed in a torus, do allow relatively 
unobscured views of the quasar core, allowing unprocessed optical flux to escape. Even 
at the most extreme viewing angles, however, current models still underestimate both the 
infrared and optical flux by several factors. 

Further modeling of the geometry and physics of the emission region in APM 
08279+5255, coupled with the effects of possible differential magnification due to 
gravitational lensing, is therefore required. Such modeling will provide a more accurate 
determination of the dust mass in APM 08279+5255, which, when coupled with 
spectroscopic observations of warm molecular gas tracers such as CO, will shed light on 
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the interaction between star-formation, warm dust and the AGN core in explaining the 
observed properties of this phenomenal system. 

It is interesting to note that pronounced infrared emission has been observed in a 
number of Seyfert Is and lis in the nearby universe. ( [Bonatto fc Pastoriza 1997|) . Similar to 
the high redshift ULIRGs, the source of this emission is thought to be a warm, dusty torus, 
with viewing-angle dependent obscuration of the continuum source and broad emission line 
region accounting for the difference between the Seyfert classes. The similarity between 



this model and that proposed for more distant and powerful ULIRGs (Parvainis et al 



1995| ), especially considering that several Seyfert I systems exhibit absorption features 



indicating the presence of high- velocity outflows [e.g. Markarian 231 flb'orster et al. I9"9"5)] 



suggests that, rather than being unique objects displaying unusual properties, systems such 
as APM 08279+5255 and H1413+117 appear to be simply more luminous members of the 
ULIRG/AGN family. 
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Fig. 1. — The SED for APM 08279+5255. The open circles present the data presented in 
this paper. The filled circles are detailed in [Irwin et al. 1998| , with the arrow indicating an 
upper-limit. The two curves are the SED for pure blackbody emitters, the solid curve for 
a system of temperature T=220K, and the dot-dashed curve for T=120K. The curves have 
been normalized to the submillimeter observations presented in this Letter. 



Fig. 2. — The submillimeter to optical SED of APM 08279+5255. The dot-dashed curve 
represents the SED for a quasar source embedded within a spherical distribution dust 
( Efstathiou fc Rowan- Robinson 1995 ), while the solid curve represents a face-on view of 
a quasar at the center of a dusty torus QGranato et al. 1996|) . The dotted line represents 
this latter model viewed in the equatorial plane. The curves have been normalized to the 
submillimeter observations presented in this Letter. Also included are the emission-frame 
SEDs for H1413+117 (triangles) and IRAS FSC 10214+4724 (crosses), normalized to the 
submillimeter SED of APM 08279+5255. As can be seen, none of the models adequately 
represent the SED of APM 08279+5255, while those of Granato describe SEDs of the other 
ultraluminous systems. 
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Table 1. Submillimeter photometry of APM 

08279+5255. The bottom line presents the 
variance-weighted means of the observations. 



Date 1350^m 850^m 450yum 



4 th April — 75±4mJy 203 ± 51mJy 

18 th April 26±4mJy 74 ± 9mJy 260 ± 130mJy 
19 th April 23 ± 3mJy 

24 ± 2m Jy 75 ± 4m Jy 211 ± 47mJy 



